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ABSTRACT 

We present a derivation of a generalized optimally-weighted estimator for the weak 
lensing magnification signal, including a calculation of errors. With this estimator, 
we present a local method for optimally estimating the local effects of magnification 
from weak gravitational lensing, using a comparison of number counts in an arbitrary 
region of space to the expected unmagnified number counts. We show that when 
equivalent lens and source samples are used, this estimator is simply related to the 
optimally-weighted correlation function estimator used in past work and vice-versa, 
but this method has the benefits that it can calculate errors with significantly less 
computational time, that it can handle overlapping lens and source samples, and that 
it can easily be extended to mass-mapping. We present a proof-of-principle test of 
this method on data from the CFHTLenS, showing that its calculated magnification 
signals agree with predictions from model fits to shear data. Finally, we investigate 
how magnification data can be used to supplement shear data in determining the 
best-fit model mass profiles for galaxy dark matter haloes. We find that at redshifts 
greater than z ~ 0.6, the inclusion of magnification can often significantly improve 
the constraints on the components of the mass profile which relate to galaxies’ local 
environments relative to shear alone, and in high-redshift, low- and medium-mass bins, 
it can have a higher signal-to-noise than the shear signal. 
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1 INTRODUCTION 

Weak gravitational lensing has proven to be a powerful 
probe of the distribution of dark matter in the Universe, 
allowing the dark mass within clusters to be mapped, pro¬ 
viding useful measurements of the halo masses of galaxies 
at a variety of redshifts, and helping to constrain funda¬ 
mental cosmological parameters, among other applications. 
However, these advances have relied only on shear, which is 
just one of the two effects weak gravitational lensing has on 
background galaxies; the other being magnification. 

When a background object is affected weakly by lens¬ 
ing, its shape will be stretched tangentially to the direction 
to the lens, known as shear, and its size will change, known 
as magnification. For measuring lensing, shear presents two 
key advantages over magnification: It has an expectation 
value of zero in the absence of lensing, and it typically has a 
larger signal-to-noise ratio than magnification. The zero un- 
lensed expectation value of shear has many benefits which 
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make measuring it easier. It allows shear to be measured 
without calibration to determine the unlensed expectation 
value, and it also indirectly leads to the result that most bi¬ 
ases in measuring shear are multiplicative rather than addi¬ 
tive. Any measurement of magnification, on the other hand, 
must deal with a comparison to a quantity that has a non¬ 
zero expectation value, which allows additive biases to play 
a role. As such, magnification requires much more care in 
analysis in order to determine an unbiased estimator. 

A relatively simple calculation of signal-to-noise which 
assumes comparable shear and convergence (as done by eg. 
ISchneider et all 1200011 can show that shear is generally ex¬ 
pected to provide much stronger signal-to-noise than mag¬ 
nification at low redshift. However, the benefit from mag¬ 
nification is not entirely negligible, and no extra observing 
time is needed to gain the information required for a mag¬ 
nification measurement; any observation which can measure 
shear can also measure magnification. As such, if an unbi¬ 
ased measurement of it can be made, it will at least provide a 
slight improvement to weak lensing mass measurements. The 
benefits of magnification are more significant at higher red- 
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shifts, due to the fact that the number of galaxies sufficiently 
resolved such that their shapes can be measured drops off 
faster with redshift than the total number of observed galax¬ 
ies. For lens galaxies at redshifts much greater than ~ 1, 
the shear signal is negligible, but the magnification signal 
is still sufficie nt to gain some information on the masses of 
these galaxies (IVan W aerbckc ct ah 20 fol: [Hildebrandt et al.l 
l2013h . Magnification can also provide additional information 
if the convergence and shear in a region of space are not of 
the same magnitude. For instance, as shear is only sensitive 
to changes in overdensity, a region of constant projected 
mass density could have positive convergence but no shear, 
allowing magnification to be used to break the “mass-sheet 
degeneracy” that prevents shear alone from determining this 
overdensity. 


It is thus well-motivated to further investigate the uses 
of weak lensing magnification in measuring the mass dis¬ 
tributions of dark matter haloes. This has been shown 


to be possible in recent works fee. Scranton et al. 2005: 

Hildebrandt et al. 2009 

: Menard et al. 2010): Ford et al. 

201U Hildebrandt et al. 

2013: Ford et al. 20141), although 


the results so far are limited, particularly compared to the 
wealth of results which have been obtained from shear data. 
These recent measurements of magnification have frequently 
used an optimally-weighted correlation function es timator 
llMenard fc Bartelmannl 120021 ; IScranton et al.l 120051) , which 
compares the projected number density of source galaxies 
in the vicinity of lens galaxies to a random catalogue, prop¬ 
erly weighted by a function of the luminosities of the source 
galaxies to optimize the magnification signalQ The mea¬ 
sured correlation function is then compared to the corre¬ 
lation function predicted from models. While this method¬ 
ology is sufficient for purposes such as galaxy-galaxy lensing, 
it cannot easily be applied to other lensing applications such 
as mass mapping. It also has the drawbacks that it requires 
significant computational time to calculate errors through 
a jackknife or bootstrap approach, and it requires a clear 
separation between lens and source samples to avoid con¬ 
tamination from clustering. 

In this paper, we start from the assumption that we 
desire an optimally-weighted estimator for the magnifica¬ 
tion itself, as opposed to previous work which has used an 
optimized correlation function estimator, and we present a 
derivation of such an estimator. Using this estimator, we 
present a new method for measuring the magnification sig¬ 
nal which avoids some of the drawbacks of the optimally- 
weighted correlation function, and we show how their result¬ 
ing estimators are related. We present a proof-of-principle 
of this new method, and we look at the results of this test to 
analyse the possible benefits magnification data may provide 
in supplementing shear data. 

In Section [2] of this paper, we present a derivation of 
an estimator for magnification and a derivation of errors 
for it, and we discuss how it can be implemented in prac¬ 
tice for galaxy-galaxy lensing. In Section [3] we discuss the 

1 There is also notable work using a conceptually-different ap¬ 
proach to magnification, relying on the observed propert ies of 
g alax ies to predict their intrinsic luminosities (eg. iHuff fe Graves! 
2011). It remains for future work to investigate if it is possible to 
combine the results of this approach with the type of approach 
presented in this paper. 


data we use for testing our method, which we take from the 
CFHTLenS. In Section [4] we demonstrate that our method 
can indeed detect magnification signals, and we look at 
the benefits of supplementing shear data with magnification 
data for galaxy-galaxy lensing. In Section [5] we discuss vari¬ 
ous aspects of the work in this paper, including a comparison 
of our method and the optimally-weighted correlation func¬ 
tion method, and we present a proposed implementation of 
our method for mass-mapping. In Section [S] we summarize 
and present our conclusions. 

In order to facilitate future work with magnification, we 
make the code used for this paper publicly-available. Details 
of its access can be found in Appendix 0 

For relevant calculations in this paper, we assume a 
cosmology with Ho = 70 km s _1 Mpc -1 , = 0.288, 
£I r = 8.6 x 10 -5 , and Uk = 0. Unless stated otherwise, all 
mass values are in units of solar masses and all distance 
values are in units of kpc. 


2 MAGNIFICATION MEASUREMENT 
METHODOLOGY 

In this section we present a derivation of a generalized 
optimally-weighted estimator for magnification and our pro¬ 
posed method for measuring it. In Section l2Tl we present the 
derivation of this estimator and show how its errors can be 
determined. In Section we present our method for mea¬ 
suring magnification for the case of galaxy-galaxy lensing. 
Finally, in Section \2 . 3 1 we discuss the key step of determin¬ 
ing the background, unmagnified source density. 

2.1 Determining an Optimal Local Estimator for 
Magnificat ion 

Let us start with the assumption that we wish to measure 
the mean magnification /r in a patch of the sky where p is 
expected to be approximately constant, such as in a circu¬ 
lar annulus around a lens galaxy. We first bin the galaxies 
observed in this patch at redshifts greater than that of the 
lens by magnitude, using small bins of width dm, which we 
will later decrease to the limit of zero. Let 

dm = n(mi)dm (1) 

be the number of galaxies observed in the i th magnitude bin 
in this patch of sky, and let 

dno,i = no(m.i)dm. (2) 

be the expected number of galaxies in this magnitude bin, 
in a patch of sky of the same area, and in the same redshift 
range in the presence of no magnification. If we assume that 
the number count can be locally approximated by a power 
law with slope 

Oii = 2.5-r— log TIO(jYli) , (3) 

dmi 

then we expect the observed number count of galaxies in 
this magnitude bin to be: 

{dm) = (4) 

llNaravanlll989l : iBroadhurst et al.lll995l) . 
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In the weak lensing limit, we expect (fi — 1) 1, and 

so we can use the approximation: 

(dm) = dno,i (1 + (u — l))“ i_1 

« dn 0 ,i (1 + (ai - 1)0 — 1)). (5) 


forms which can be constrained by observations, also giving 
us a functional form for w(m). For a set of observed galaxies 
with magnitudes we can express the differen¬ 

tial number count n(m) in terms of a sum of Dirac delta 
functions: 


Solving this equation for ji, we get: 


ji — 1 -t 


(dm) — dno,i 

(ai l)d?To,i 


( 6 ) 


If we assume the statistical error on estimating dno,i and a; 
to be significantly less than the standard error of dm, we 
can then use 


m — 1 + 


dm — dno.i 
(at — 1 )dn 0:i 


( 7 ) 


as an estimator for fi, with (fn) = U, and which is indepen¬ 
dent for each magnitude bin i. 

In order to optimally combine the estimators fn for all 
magnitude bins, we apply inverse-variance weighting. We 
start with the standard deviation of each estimator, which, 
under the assumption that dn; is Poisson, can be calculated 
to be: 


~ , 1 1 - • (8) 

| a t - 1 | ^/dn 0 ,i 

The derivation of this equation is presented in Equation (EH) 
in Appendix m This then gives us the optimal inverse- 
variance weighting for bin v. 


dlVi 1 opt — 



dno,i (at — l) 2 

1 + (ai - 1 ) (fi — 1) 


( 9 ) 


While the weighting scheme in Equation © will provide 
the estimate of fi with the minimum standard error, it re¬ 
quires foreknowledge of the actual value of fi to determine 
the weights. For a likelihood-based approach, we can use 
the value of fi predicted by the model for the weight calcu¬ 
lation, but this will require a separate summation over bins 
to calculate the likelihood for each model ji. Since we expect 
(/.i — 1) <C 1, we instead use the weighting scheme: 


dwi = dno,i (at 


l) 2 


( 10 ) 


which is model-independent. This will save significant com¬ 
putational time at the expense of a negligible decrease in 
the precision (but not accuracy) of our constraints. 

We can then calculate our estimate of /.i for this patch 
of sky through: 


, _ E idwi/n 

M E; dwi 

This will have standard error: 



( 11 ) 


( 12 ) 


the accuracy of which depends on the accuracy of the stan¬ 
dard deviation we calculated for each m in Equation ©. 

Let us now examine the case where the bin size dm ap¬ 
proaches zero. We now express the parameters drii, drio.i, 
ai, fu, and dwi as dn(m), dno(m), a(m), /x(m), and dw(m) 
respectively. Here, m is the magnitude of the differential 
bin, and we define dn(m) = n(m)dm, dno(m) = no(m)dm, 
and dw(m ) = w(m)dm. no(m) and a(m) have functional 


n(m) = E 5(m — mj). (13) 

3 


This form gives the expected behaviour that integrating 
n(m) from m to m + dm results in the number of galax¬ 
ies with magnitudes in this range. 

Let a represent the minimum observable magnitude for 
galaxies, and let b represent the maximum magnitude used 
in the catalogue. We can then calculate our estimate of fi 
through: 


wL 


= i + 


w(m)[i(m)dm. 

1 f b 

— / \n(m) — no(m)] [a(m) — 1] dm (14) 

v J a 


where: 


rb rb 

W= w(m)dm = / no(m) (a(m) — l) 2 dm. (15) 

J a J a 

Expanding n(m) using Equation (| 1 3 1) . this gives us: 

jl = 1 + 

(E H m 3') - !] - J n o( m ) [a(m) - 1] ( 


— 1] dm. , (16) 


where mj represent the magnitudes of observed galaxies in 
the range a < mj < b. 

As before, this estimator will have the standard error: 

ag « W~ 1/2 , (17) 


and will be approximately Gaussian for large n. The ac¬ 
curacy of this depends on the validity of our assumptions 
on the distribution of number counts used in Equation ©. 
In practice, an empirical error estimate (determined, for in¬ 
stance, from the standard deviation in the fi for various sim¬ 
ilar regions) is likely to be more reliable. 

From this, we can estimate the likelihood for a given 
value of fi to be 


C(ji\fi) 



: ^i J 

W(u-u) ^ 


(18) 


The weighting used for this estimator is in fact the 
same as that generally used for the optimally-weighted 
correlation function, and w as originally presented by 
iMenard fe Bartelmannl (l2002i 'l without derivation. We will 
later show in Section O that this estimator is simply re¬ 
lated to that used for correlation-function-based methods 
in the scenario in which the lens and source samples are 
equivalent. In the following section, we will show how the 
estimator presented here can be used as a basis for an al¬ 
ternative method for calculating a magnification signal by 
directly comparing the observed source counts to the mea¬ 
sured background density at a given redshift. 
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2.2 Implementation for Galaxy-Galaxy Lensing 

We start with the case of a galaxy-galaxy lensing analy¬ 
sis, which will allow us to readily test our implementation 
against the lensing signal measured with shear information. 
We start by binning our lens sample by both stellar mass 
and redshift so that the lenses in each bin will have approx¬ 
imately the same halo mass profiles. 

For each lens, we bin the possible patch of sky in which 
sources could potentially appear based on projected distance 
at the redshift of the lens. This gives us a series of annuli 
for each lens, some of which may overlap masked regions 
of the image. We determine the unmasked fraction of each 
annulus through the ratio of unmasked pixel centres to all 
pixel centres inside the annulus. In practice, most masks 
are relatively coarse, and so this fraction can be reasonably 
estimated by regularly sampling pixels in order to conserve 
computational time. 

The angular area of the k th annulus around the Z th lens 
is then: 


flfei = A 2 /fci7T (R 2 + 1 — Rl) , (19) 

where A is the angular diameter distance at the redshift of 
the lens, fki is the unmasked fraction of this annulus, and 
Rk and Rk+i are its inner and outer radii respectively in 
physical units. 

To ensure minimal contamination of the source sample 
with galaxies which are actually in the lens plane, we require 
any sources used for our analyses to be separated from their 
associated lenses by a redshift buffer A z, which we set to 
0.2. As galaxies in the lens plane are clustered with each 
other, any which are misidentified as sources will be clus¬ 
tered with lenses as well, resulting in an artificial increase 
in the observed counts of sources around lenses, biasing the 
magnification signal. It is therefore necessary that this buffer 
be sufficiently large that this contamination has negligible 
effect on the measured magnification signal. See Section [3T] 
for discussion of our choice of Az. 

Note that this separation between lenses and sources is 
applied to each lens individually. This allows for the pos¬ 
sibility that the lens sample might overlap with the source 
sample. Previous methods for measuring magnification have 
relied on the two samples being well-separated. The loos¬ 
ening of this requirement for our method will allow for an 
increased signal-to-noise ratio in situations in which it is 
beneficial to use overlapping lens and source samples. 

The expected unmagnified number count of sufficiently- 
distant sources in a given magnitude bin projected in this 
annulus is 


no,ki(m)dm — nun(m, z + Az)dm, (20) 

where 2 is the redshift of the lens and no(m, 2 + A z)dm is 
the number of sources per unit area with magnitude between 
m and m + dm at redshift z + Az or greater. This gives us 
the relevant no(m) to use in Equation and a(m) can 

be determined through Equation ©, allowing us to get the 
estimate Afci for this annulus. 

This process can be performed more efficiently through 
storing the integrals 

d(z) = r n(m, z) [a(m) — 1] dm (21) 


and 


rb 

w(z) = / n(m, z) [a(m) — l] 2 dm (22) 


for all needed redshift values and calculating Afci through: 


Afci = 

1 + IWzl + A2) (E [“to”) - !] - + A ^) ’ 

(23) 

where mjki represent the magnitudes of all sufficiently- 
distant sources observed in this annulus, and a and b rep¬ 
resent the lower and upper magnitude limits of the source 
sample, respectively. 

By comparing Equation m to Equation m and 
Equation <tm we can estimate the standard deviation of 
Afci by 

o-j i fe! = [fluw(zi + A2)] _1/2 . (24) 

When we stack lenses, we can use this estimate to apply 
inverse-variance weighting, giving each individual annulus a 
weight of 

Wki = 0.kiw(zi + Az). (25) 


Finally, we get the estimated magnification for the 
stacked lenses in this annulus from the weighted mean across 
all lenses: 


j^k 


Ei w kldkl 
Ei W kl 


(26) 


If our estimates of standard deviations in number counts 
were perfectly accurate, the standard error in this estimator 
would be 


Optfc ,ideal 



(27) 


but as this is unlikely to be the case (for instance, due to the 
clustering of sources making their distribution not perfectly 
Poisson), we instead determine it empirically. We start by 
empirically determining the weighted standard deviation in 
the jlki estimates for the set of lenses, 


k — Ck 


Ei WkiPh A E; Wkifiki Y 

T,i w ki V E i w ki ) 


where 


(28) 


C k = . J2lWkl = (29) 

VlEiWfci ) 2 -Ei^ii 

rtThePawnlbOliM We then estimate the standard error in 
the magnification estimate Afc f° r this set of annuli through: 


cr — 


a P-k 

VNk-T 


(30) 


where Nk is the number of lenses with non-zero weight for 
this annulus. 


2 This formula for the standard error of a weighted mean was 
chosen based on a comparison of various estimators, as it was 
found to generally have the lowest amount of bias in various test 
cases, some of which involved biased weights being used. 
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This estimate of the magnification can then be con¬ 
verted into an estimate of the convergence: 

= 1 - s/hk 1 +7t,fci ( 31 ) 

where 7t,fc is the tangential-shear estimate for this annulufl 
The standard error in this estimate is 


&Pi,k - 


07i 2 + 27t,feCT^2 


(32) 


^ 4 [fik 1 + 7t,fc) 

The estimated projected overdensity in this annulus is then 
Efc — E C rit ,AAfc, (33) 


where E cr ; t ,fc is the critical density for this annulus, 


E 


crit,fc 


c 2 A 
47tG c A As ’ 


(34) 


where A is the angular diameter distance from redshift zero 
to the mean source redshift, A from redshift zero to the 
mean lens redshift, and As from the mean lens redshift to 
the mean source redshift. The standard error in the pro¬ 
jected overdensity is then: 


a 




S 


crit,fc^&,fc • 


(35) 


2.3 Determining Umnagnified Source Density 

In measuring the magnification signal of a given region of 
space, we require knowledge of what the source density 
would be as a function of magnitude in the absence of any 
magnification. In Equation O and Equation dm this rep¬ 
resents no(ra), which is used to determine a(m) through 
Equation J3]) and W through Equation (1151) . Given the large 
area covered by the CFHTLenS, it is a relatively safe as¬ 
sumption that the survey region is on average unmagnified. 
There will be some overdense and underdense regions, but 
the region as a whole should have a projected density ap¬ 
proximating that of the whole sky, and so we can thus use 
the measured source densities across the entire field to de¬ 
termine no(m). 

This approach has two key advantages over a theoreti¬ 
cal model. First, it accounts for many observational effects 
which will identically affect the source samples we later use 
for measuring magnification, such as the difficulty in identi¬ 
fying faint galaxies. Second, as we are interested in all galax¬ 
ies which appear behind a given annulus, we are in fact inter¬ 
ested in the number density of all galaxies at a given redshift 
or greater, whereas models are typically fit only to galaxies 
within a single redshift slice. This resulted in us being un¬ 
able to satisfactorily fit a model such as a Schechter function 
to the observed galaxy counts beyond a given redshift, even 
allowing for perturbations to the functional form. 

This approach does result in noise being present in 
no(m). As it is measured over a much greater area of the sky 
than even the stacked regions we wish to measure magnifi¬ 
cation for, however, this noise is likely to be sub-dominant 

3 Equation ED uses the assumption that the cross-shear 7 X ^ is 
expected to be zero, as is the case in galaxy-galaxy lensing. For 
applications where this is not the case, 7 /- = y 7 ^ + 7 ^ k should 
be used in place of 7 t & in this equation and Equation 1321) . 


to the Poisson noise in n(m), the number count distribution 
of galaxies which actually appear in these regions. However, 
we must differentiate no{m) to determine a(m), and differ¬ 
entiation amplifies this noise. In order to handle this, we first 
smoot h the measured no(m) using the ISavitzkv fe Golavl 
( 4963 ) algorithm prior to differentiation. 

As can be seen most easily in Equation da, the esti¬ 
mator for magnification is particularly sensitive to no- For 
instance, in the simplified scenario where a(m) = 2, a bias 
of 0.1 per cent in no leads directly to a bias of of —0.001 
in (i. As the typical magnification signal is of this order, we 
must take care to minimize bias in no as much as possible. 

If the positions of lenses and sources were completely 
independent, the approach we use would be expected to be 
unbiased when applied to galaxy-galaxy lensing. However, 
even aside from contamination of the source sample with 
galaxies in the lens plane misidentified as sources, some cor¬ 
relation or anti-correlation between the samples is expected. 
This can be caused by a number of position-dependent ef¬ 
fects which can introduce biases into the redshift estimates 
for galaxies, including atmospheric seeing, depth variations 
due to the dithering strategy of the survey, and depth vari¬ 
ations due to galactic dust. 

The method use d for estimating photome tric redshifts 
in the CFHTLenS bv lHildebrandt et al.l J2012I ') used a local 
model for atmospheric seeing to minimize the bias in galaxy 
colours due to it, so this is unlikely to be an issue. However, 
regions with poor seeing and regions with poor depth will 
have larger errors on the estimated colours of galaxies. If 
there is a gradient in the object density in colour space (eg. 
more blue galaxies than red), then noise in colour determi¬ 
nations will result in asymmetric scatter along this gradient. 
Even if any particular galaxy is as likely to be scattered to 
a higher colour as a lower colour, the net result will be to 
smooth out the object density distribution in each axis of 
colour space. This will then result in the estimated redshifts 
of these objects being systematically biased. As the relation 
between redshift and colour is not linear, the direction and 
magnitude of this bias will vary with redshift. This means 
that some redshift bins may be spuriously underdense, while 
others are spuriously overdense. 

Let us present an example of what might occur, and 
how this might lead to a systematic bias in the magnifica¬ 
tion signal. Consider an extreme example: a field is divided 
into two regions of equal size; one region has good seeing, 
and another region has poor seeing. This poor seeing re¬ 
sults in the measured luminosities of galaxies in this region 
having greater errors than the luminosities of galaxies in the 
other region. If we suppose that there are more blue galaxies 
than red and that, for the sake of this example, the distri¬ 
bution between the two extremes is smooth, then we would 
expect more galaxies to be scattered to redder colours than 
we would expect to see galaxies scattered to bluer colours. 
Thus, there would be a net bias in this region for all galaxies 
to appear redder. Let us assume a simple photometric red¬ 
shift estimation which relies only on this colour. In this case, 
this would result in a systematic bias in this region for all 
galaxies to have higher redshifts. As there is more volume 
at higher redshifts, there are more galaxies there. Thus, a 
bias of all galaxies to higher redshifts would result in a lower 
density of galaxies in any given redshift bin. 

We thus have one region with good seeing and approx- 


© 2015 RAS, MNRAS Q00.mi20l 
















6 B.R. Gillis and A.N. Taylor 


03 

C5 



5000 6000 7000 8000 9000 

Galaxies per deg 2 in 0.3 < z < 0.4 


10000 11000 


Figure 1 . A comparison of the number densities of galaxies in 
two different photometric redshift slices in the CFHTLenS. Each 
point represents an individual field. The anti-correlation between 
the two redshift slices implies the presence of a systematic error 
in redshift estimates which varies from field to field. 


imately normal density of lenses and sources, and one with 
poor seeing and a lower density of both lenses and sources. In 
measuring no, we would weight both of these regions equally 
and calculate an expected source density which is the mean 
of the two. However, in measuring the magnification signal, 
we would be using the lower number of lenses in the region 
which also has a lower number of sources. Our measured n 
would thus be more akin to a weighted mean, which in this 
case would end up being weighted more by the region with 
a lower density of sources. Thus, even without any magnifi¬ 
cation, we have a scenario in which the measured count of 
galaxies differs from the expected count. 

Evidence of this type of bias can be seen in Fig. [T] which 
shows the galaxy number densities of different fields from the 
CFHTLenS in the redshift slices 0.3 — 0.4 and 0.8 — 0.9 plot¬ 
ted against each other. This plot shows clearly that fields 
which have greater densities in the 0.3 — 0.4 slice have lower 
densities in the 0.8 — 0.9 slice, and vice-versa. There are 
also many case of well-separated slices which have corre¬ 
lated number densities, but the anti-correlation here is most 
striking, as it is difficult to explain with any effect other 
than non-trivial redshift systematic errors such as we have 
discussed. More insight into thi s can be gained from th e 
rightmost panel of Figure 4 of iHildebrandt et ah] (|2012l h 
which shows a comparison of the photometric redshifts as 
used here against the spectroscopic redshifts of a matched 
sample of galaxies, and observing the complicated shape of 
the distribution. The field-by-field comparison we show here 
demonstrates that this distribution is not constant between 
fields. 

A full discussion of this type of effect is beyond the scope 
of this paper. Interested readers are directed to Morrison and 
Hildebrandt (2015, submitted), which discusses this issue in 
detail. 

In the example above of how this effect might lead to 
an error in the measured magnification signal, we pointed 
out that the result is similar to the difference between the 
unweighted and weighted means of source densities. We can 


thus correct for it by weighting our initial calculation of 
mean source density by the density of lenses in the same 
region of space. The more local this weight is, the noisier 
it will be, and the more it will become entangled with the 
actual magnification signal, but the better it will account 
for this effect. 

For this paper, we have decided to weight on a field-by- 
field basis, as the discontinuity in both position and time of 
observation between fields is likely to result in much more 
variation in seeing and depth than the variation across the 
field-of-view within a single field. To do this, when we are 
counting all galaxies in an individual field at a redshift z B or 
greater, we weight this by the number of lenses at redshif'tQ 
z s — A z. See Section EU for discussion of the impact of this 
weighting on our results. 


3 DATA 

In order to test our methodology, we use the 
publicly-available galaxy catalogues from the 
Canada-France -Hawaii Telescope Lensing Survey 
(iHevmans et al.1 l2012h . hereafter referred to as 
“CFHTLenS,” which are available for download from 
http://www.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/ 
community/CFHTLens/query .html. The CFHTLenS is 
a 154 deg 2 survey (125 deg 2 after masking) llErben et ali 
l2013h . based on the Wide component of the Canada-France- 
Hawaii Telescope Legacy Survey, which was observed in 
the period from March 22nd, 2003 to November 1st, 2 008, 
using the MegaCam instrument (iBoulade et al.l l2003l h It 
consists of deep, sub-arcsecond, optical data in the five 
optical SDSS-like filters u* g'r'i' z'. The data is divided into 
171 individual patches of sky, which we refer to as “fields” 
in this paper. 

In order to avoid contamination from foreground stars 
and galaxies, cosmic rays, and other observational de¬ 
fects, masks were applied to the CFHTLenS fields usin g 
the automask tool (iDietrich et al.ll2007l : lErben et alll2009h . 
Shapes of background galaxies with i' < 24.7 in the un¬ 
masked regions were measured with the lens fit shape mea¬ 
surement algorithm (iMiller et alj|20l3 l. giving an effective 
density for sources with sha pes of 11/arcmin 2 in the red¬ 
shift range 0.2 < z p o ot < 1.3 llHevmans et al.ll2012lh For the 
analysis in this paper, we use all fields in the survey, not sim¬ 
ply those that passed the systematic s tests for cosmic shear 
measurements (jHevmans et al.ll2012i i. as it has been demon¬ 
strated that fields with systematics that may affect cosmic 
shear have no effect o n galaxy-galaxy lensing measurements 
(IVelander et al.|[2014l L 

Photometric redshi fts for the survey were estimated 
with the BPZ code jBemtedl2000l h The accuracy of photo-z 
measurements for the CFHTLS-Wide survey was improved 
through homogenizing the PSF through different bands in 
the CFHTLenS survey. Photo-zs were made available for 
the entire survey, with a typical redshift uncertainty of 
~ 0.04(1 + z) in the r edshift range 0.2 < z p hot < 1-3 
llHildebrandt et ahll2012l L 


4 The photometric redshifts we use are binned to discrete inter¬ 
vals of 0.01. hence the precise value being used here. 
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We use the st ellar mass estimates described by 
IVelander et al.l (I2014I ). obtained by fitting spectral energy 
distributi on (SEP ) templates, following the method of 
illbert et al.l (l2010l ) . These stellar masses were found to be in 
general agreement w ith deeper data su ch as WIRDS, which 
includes NIR filters llBielbv et al.ll201~3) . up to z = 0.8. 

We use the magnitudes measured with the r' filter for 
all galaxies. We considered using the i'-band magnitudes, 
but this was found to be problematic, as the number den¬ 
sity of galaxies as a function of i'-band magnitude in the 
CFHTLenS shows a prominent step at mp = 23, which 
would result in a(m), the logarithmic derivative of this func¬ 
tion, being discontinuous at this value. This is possibly an 
indirect result of the fact that the original i' filter used for 
the CFHTLenS had to be replaced partway through obser¬ 
vations. The r'-band magnitudes showed no such issue, and 
typically had the next best signal-to-noise after the i'-band 
magnitudes, and so we decided to use them. From here on, 
whenever we refer to magnitudes, we will thus be referring 
specifically to the r'-band magnitudes. 


3.1 Sample Selection 

We select galaxies from the public catalogue using the fol¬ 
lowing cuts on it: 

• Z_B < 4 

• MAG_R < 24.7 

• starjlag < 0.01 

• CHI_SQUAREDJBPZ < 2. 


This provides us with a sample of galaxies with reasonable- 
quality redshifts. Other cuts were tested to improve redshift 
quality, such as imposing a cut on the “ODDS” parameter, 
but they were found to cause too much bias in the resulting 
galaxy selection. We use Z_B as the photometric redshift 
estimate for each galaxy. 

We do not impose a cut on the “MASK” parameter. 
Rather, we download the mask files and apply them our¬ 
selves, removing any galaxy from the sample which resides 
in a region with MASK > 1. We do this as we found that 
the catalogue’s MASK values for a small number of galax¬ 
ies did not agree with our determinations, possibly due to a 
different version of the mask files having been used to assign 
the values in the catalogue. Assigning the MASK values our¬ 
selves thus ensures consistency with the mask files we have 
available, which we also use for determining the unmasked 
fractions of annuli around lenses. 

We use the entirety of these galaxies for our source sam¬ 
ple, and all galaxies with Z_B < 1.3 for our lens sample. In 
selecting lens-source pairs, we require that the source be in 
the background of the lens, separated by a minimum red¬ 
shift of Az = 0.2. This value is chosen to minimize the the 
fraction of galaxies identified as appropriate sources which 
in fact reside in the same bin as the lenses. To determine 
this fraction, we matched the CFHTLenS galaxy catalogue 
to publicly-available spectroscopic galaxy catal ogues f rom 
the VIMOS VLT Deep Survey (“VVDS” iLe Fevre et all 
120051; Garilli et al.l 200 8li and th e DEEP2 galaxy survey 
( Davis et al .. 2007T iNewman et al . 20131) , both of which over¬ 


lap portions of the CFHTLenS. We then calculated the frac¬ 
tions of galaxies with sufHciently-large photometric redshifts 
which have spectroscopic redshifts within a given bin. As 



Lens Redshift Bin Centre 


Figure 3. For a given lens redshift bin of width 0.1, the fraction 
of galaxies with photometric redshifts at least a given redshift A 2 
greater than the centre of the bin (corresponding to the source 
sample), which are likely to actually reside in the bin. We deter¬ 
mine this fraction through matching galaxies in our sample with 
spectroscopic redshift catalogues from the VVDS and DEEP2, 
and determining the fraction of galaxies with sufficiently-large 
photometric redshifts which have spectroscopic redshifts within 
a given bin. Larger Az values consistently correspond to smaller 
contamination fractions, though this effect seems to converge to 
a maximum by Az ~ 0.2. The contamination fraction values for 
each lens redshift bin are plotted near the centres of the bins, 
with a horizontal offset applied for visual clarity. 


shown in Fig. [3] a buffer of Az = 0.2 ensures that fewer 
than 1 per cent of galaxies identified as appropriate sources 
in fact reside in the same bin as the lenses for lens redshift 
bins Zmid < 0.65, and the benefit of increasing this value 
further is indistinguishable from noise. 


4 APPLICATION TO DATA 

In order to provide a proof-of-principle that our proposed 
method can be successfully used to measure the magni¬ 
fication signal from weak gravitational lensing, we test 
it by measuring the galaxy-galaxy lensing signal in the 
CFHTLenS. We choose the galaxy-galaxy lensing applica¬ 
tion as it is computationally straightforward to measure the 
shear signal at the same time as the magnification signal. 
As weak lensing shear is a well-studied and robust method 
for determining halo mass profiles with signal-to-noise of the 
same order of magnitude as weak lensing magnification, it 
provides an ideal comparison. 

One drawback in using shear to model halo mass pro¬ 
files is that the shear signal is degenerate to the addition or 
subtraction of a sheet of constant projected mass density to 
the held of view. Magnification is not subject to this degen¬ 
eracy, which leads to the result that a comparison with shear 
alone cannot be used to determine whether a hat offset in 
the magnihcation signal is due to a bias in the measurement 
or to the presence of an actual over- or underdensity in the 
region of interest. See the discussion on the potential for 
lens-source correlations or anti-correlations to cause such a 
bias in the magnihcation signal in Section \2 . 3 1 
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Figure 2. Number density of galaxies at a given redshift or greater as a function of r'-band magnitude (left) for the CFHTLenS, and 
the corresponding calculated values of a (right). The strong oscillations observed in a for high redshifts and low magnitudes are due to 
Poisson noise, due to the low numbers of galaxies in these bins. This results in these bins being correspondingly down-weighted, and so 
these oscillations are likely to have no appreciable effects on our results. In the left panel, the lowest redshifts correspond to the highest 
lines, and in the right panel, the lowest redshifts typically correspond to the lowest lines. 


In Section f4. 11 we discuss the model we use for the ex¬ 
pected form of the lensing signal, and in Section 14.21 we 
present the results of our tests. 


4.1 Modelling Lensing Signal 


To fit our measured lensing signals, we assume that the 
underlying mass distribution can be described by a profile 
which is on average radial, and which can be described with 
relatively few parameters. The most common mode l used 
for this purpose is know n as the halo model dSeliakl l2000l : 
iMandelbaum et al.l 12003 ). which has two free parameters: 
the halo mass of the galaxies in the sample, and the satel¬ 
lite fraction, with other parameters directly fixed through a 
physical model. For this paper, however, we choose to use 
a different model, which has one additional parameter to 
characterize the surrounding environment of galaxies in the 
sample. This allows us to better investigate how magnifi¬ 
cation might aid in breaking degeneracies between separate 
factors which characterize the environment of galaxies. 

Our model uses four free parameters: Mih, the total 
one-halo mass of the lens; M gr , the mean mass of the group 
or cluster a typical satellite lens resides in; / aa t, the frac¬ 
tion of lenses which are satellites of a group; and Koffaet, a 
flat offset to the measured k due to observational biases. 
Other relevant parameters are fixed directly from observa¬ 
tions or with fitting functions. We model the galaxies in 
each bin as having thei r mass domina t ed b y a dark matter 
halo with a truncated Navarro et al.l dl997f) profile, using 
the model of lBaltz et alJ ( 20091 ) (hereafter a “tNFW” pro¬ 
file). We leave the mass of this profile as a free parameter, 
determine t he con centration through the fitting function of 
iNeto et al.l J2007I ). and set the truncation radius equal to 
double the virial radius. 

Of these galaxies, we assume that a fraction / sa t are 
satellites of a group, where / sa t and M gr , the total mass of 
the group, are free parameters. We assume the group’s total 


mass can also be modelled with a tNFW profile, and we de¬ 
termine its concentration and truncation radius in the same 
manner as for the galaxy’s halo. We assume that the distri¬ 
bution of satellites relative to group centre can be described 
by a scaled tNFW density profile with concentration c = 2.5 , 
which is consistent with the results of lMcGee et alJ l|2009l ). 
The details of calculating the contribution of t he group en¬ 
vironm ent to the lensing signal can be found in lGillis et al.l 
d2013h . which uses the same model for shear except with a 
fixed /sat. For computational simplicity, we do not include 
any contributions from large-scale structure to the lensing 
signal (the so-called “two-halo term”) at this stage of anal¬ 
ysis, but we plan to do so in future work. 

This model results in a mean mass profile for a given 
set of parameters Mih, M gr , / aa t. To determine the shear 
and magnification lensingsignals, we project this mass pro¬ 
file along the line of sighl[j to get the projected overdensity 
E (^Rj . The expected shear and convergence signals can then 
be obtained through: 


and 


(7t> = 


E (< R) - E (R) 


(k) = 


s(fl) 

Scrit 


(36) 

(37) 


where E ( R ) is the mean projected density at a distance 
R from the galaxy, E (< R) is the mean projected density 
contained within a distance R of the galaxy, and E cr it is 
the critical projected density, using the same definition as 
Equation (13411 . 

Finally, we assume that the correlation and anti¬ 
correlation between the number counts of lenses and sources 


5 For computational simplicity, we neglect the presence of struc¬ 
ture correlated to the galaxy and group, the so-called “two-halo” 
term in the lensing signal. 
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across the field of view due to observational effects will have 
the net effect of shifting the measured magnification signal 
by a constant value, which we leave here as a free parameter, 
^offset- While the other three free parameters affect both the 
measured shear and magnification signals, Koffset only has an 
effect on the magnification signal. 

We determine the best-fit set of parameters and their as¬ 
sociated errors for each bin of interest through a three-phase 
Markov chain Monte Carlo (MCMC) algorithm. In the first 
and se cond phases, we follow the Metropolis-Hastings algo¬ 
rithm ([Metropolis et al.lll953l f. using a y 2 comparison of the 
model lensing signal with the measured values to determine 
the acceptance ratio. The first phase is a “burn-in” period, 
in which the algorithm is given time to reach the region with 
the lowest \ 2 - In the second phase, we periodically record 
the points arrived at by the algorithm. The standard devi¬ 
ations in these values represent the lcr errors in the best-fit 
model parameters. In the final phase, we apply simulated 
annealing, successively decreasing the size of steps and the 
acceptance ratio, to arrive at the best-fit set of parameters. 

For each bin of interest, we fit the data in three man¬ 
ners: To magnification information only, to shear informa¬ 
tion only, and to both magnification and shear information. 

4.2 Resulting Measurements of Magnification 
Signal 

In order to test that our model is performing as expected, we 
first perform a measurement of tangential shear in various 
lens stellar mass and redshift bins and test fitting out model 
to the resulting signal. The results of this test are shown 
in Fig. [4] where we present the shear signal in terms of the 
reduced tangential shear 7 t against angular separation for 
each lens bin as well as the predicted shear for the best-fit 
model. In each bin, the model passes a \ 2 test in comparison 
with the data, as does the mean y 2 of 17.6 for 15 degrees of 
freedom and 15 bins. Our model thus appears to provide a 
good fit to the shear data. 

In Fig. 0 we show the measured magnification signal 
for the same lens stellar mass and redshift bins as in Fig. [4] 
presented in terms of the convergence ft, which is related 
to the measured magnification through Equation m- Also 
plotted on it are curves representing the predicted ft from 
the best-fit model which uses only magnification data, the 
best-fit model which uses only shear data (for which 7 t for 
each bin is plotted in Fig. 11 , and the best-fit to the com¬ 
bination of the shear and magnification datasets. As the fit 
using shear data alone does not constrain ft 0 ff se t, the fitted 
value from using magnification alone is used for that model. 
The data generally shows the expected form and is well-fit 
by our model, though the errors are typically larger than 
for the shear measurements. Notice that in some bins the 
signal converges to a positive value at large angular separa¬ 
tion between lenses and sources, and in others it converges 
to a negative value. This is most likely an effect of variable 
seeing and other position-dependent effects leading to corre¬ 
lation or anti-correlation between lenses and sources, as we 
discussed in Section and which is the reasoning behind 
our use of the ftoffset model term. 

We note that in many cases, the innermost radial bins 
show very high magnification signals compared to the mod¬ 
els. We believe this to be most likely due to the obscuration 


of background galaxies by dust in the regions of lens galaxies 
(particularly those which reside within groups or clusters). 
Dust can most easily obscure the faintest magnitude galax¬ 
ies, and it is at the faintest magnitudes that the slope of the 
luminosity function is shallowest (as can be seen in Fig. H>. 
At these magnitudes, a(m) — 1 is negative, and so the the 
obscuration of a galaxy would result in a spurious increase 
in the magnification signal. As we know that the innermost 
radial bins are likely to be problematic in this manner, we 
leave them out of our fitting procedure and calculation of \ 2 
values. Methods have been proposed to handle dust extinc¬ 
tion (eg. I Menard et al .1120 id . which finds that magnification 
is typically a stronger effect than dust extinctio n by a factor 
of ~ 3 in the V band, and iBauer et alj|20l3 l. but as this 
paper is intended to serve as a first demonstration for the 
magnification methodology proposed here, we leave the full 
application of these for future study. 

In the low redshift bins, the model which is best fit to 
only shear data generally also does a good job at fitting 
the magnification data. It fails a \ 2 t es t at the 5 per cent 
threshold for only five of the fifteen low-redshift bins, and 
the fit to the combined dataset fails for only two of these 
bins: the 9 < log Miens /< 10 and 0.3 < 2 < 0.4 bin, and 
the 11 < logMiens/Af® < 12 and 0.6 < z < 0.7 bin. This 
failure ratio of 13 per cent for the combined fit does imply 
that something beyond statistical error is to blame here, but 
this is not surprising given the simple model we use in this 
paper. 

In the higher redshift bins, the best-fit model to shear- 
only does a much poorer job at fitting the magnification sig¬ 
nal. This is not surprising, as at these redshifts the signal- 
to-noise of the shear measurements is much lower, and so 
noise in the shear signal can propagate through into noise 
into the best-fit measurements, resulting in a high \ 2 value 
when it is compared with the magnification data. The fit to 
the combined dataset does a better job at fitting to both, 
though it still fails at the 5 per cent threshold for seven of 
the eighteen high-redshift bins. It remains for future inves¬ 
tigation to determine whether this is due to issues with the 
shear signal, the magnification signal, the modelling process, 
or a combination of these factors. 

Overall, the strong fit of our model implies that our 
method is successfully recovering a real magnification signal. 
This magnification signal may, however, still have biases in 
it due to effects we haven’t fully corrected for, such as dust 
obscuration or contamination of the source sample. Thus at 
this stage, the specific values fit by our model may still be 
biased. Given the strong agreement with the shear-signal, 
these biases are likely at a small level, however, and we can 
still gain insight into how magnification data can be used to 
improve the errors on fits to various parameters. 

This can be seen in Fig. [6j where we show the best- 
fit parameters for our model fits and their associated errors 
for all lens stellar mass and redshift bins tested. This fig¬ 
ure shows the fits to shear data alone, magnification data 
alone, and to the combined dataset. As shown in Fig. [3j 
the contamination of the source sample with galaxies in the 
lens plane which were spuriously misidentified as sources be¬ 
comes a significant issue at redshifts greater than ~ 0.7, and 
so the fits including magnification data at these redshifts 
may be less reliable. The fitted M\h and M gr parameters 
from magnification alone seem to be in general agreement 
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Figure 4. Measured shear signal 7 (points and errorbars) for different lens redshift and mass bins. Also shown is the best-fit model 
using shear data alone (dashed line). The top panel displays low-redshift bins, and the bottom panel displays high-redshift bins. The y 2 
values represent the y 2 for a comparison of the shear data with the shear-only fit (15 dof.). 


with those from shear alone. It is in the / aa t parameter that 
the two fits differ, with the magnification-only fit regularly 
favouring a larger / sat . 

Interestingly, some of the fitted /toffset values show a sig¬ 
nificant change between the magnification-only fit and the 
fit to the combined dataset. This seems to occur when the 
best-fit group mass parameter M gr for magnification alone 
was spuriously high but was corrected to a more-reasonable 
value when complemented with the shear data. There are 
also bins in which the shear-only fit appears to show a spu¬ 
riously high or low, but in which the fit to the combined 


dataset is much more reasonable. This clearly shows the po¬ 
tential benefit of combining shear and magnification data. 

The trend observed in the best-fit values for the one- 
halo mass parameter Mih at large redshift is potentially 
worrisome - the combined fit seems to consistently favour 
very low masses, often lower than either the shear or magni¬ 
fication data alone favours. This is possibly due to the fact 
that we are not fitting to the innermost radial bins, which 
means we are losing a significant amount of information rel¬ 
evant to the one-halo mass fit. This could also potentially 
be due to contamination of the source sample in high lens- 
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Projected Separation (kpc) 
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Figure 5. Measured magnification signal k (points and errorbars) for different lens redshift and mass bins. Also shown are the best-fit 
models using shear data alone (blue dashed line) using magnification data alone (red solid line), and for the combination of the shear 
and magnification data (black dot-dashed line). The top panel displays low-redshift bins, and the bottom panel displays high-redshift 
bins. As the fit using shear data alone does not constrain K 0 fj se t; the fitted value from using magnification alone is used for that model. 
The x 2 values represent the \ 2 for a comparison of the magnification data with the magnification-only fit (14 dof.), the shear-only fit 
(17 dof.), and the combined fit (14 dof.), respectively. 


redshift bins leading to spuriously low magnification signals. 
We discuss this possibility further in Section 15. 3 1 

In Fig. [7] we show the ratios of the calculated errors on 
our fitted parameters to the errors obtained when shear in¬ 
formation alone is used. Both the ratios of the magnification- 
only fit and the combined dataset fit are shown, as well as 
a calculation of the theoretical ratio, where the theoretical 


combined error is: 

Uth = (f^shear + ^mag) V2 , (38) 

which is not necessarily valid for a fit of more than one 
parameter, depending on the relative degeneracies of the 
two fits. 

This plot shows that when used alone, magnification 
typically does a poorer job at constraining the model pa- 
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Figure 6. Best-fit values and errors of our model parameters, for various lens redshift and mass bins. From left to right, the columns 
show the parameters: One-halo mass, Group mass, Satellite fraction, and Ac 0 ffset- The best-fit values using shear data alone are shown 
in blue with solid circle markers, values from magnification data alone are shown in red with open circle markers, and values from the 
combined data are shown in black with star markers. As shear data alone cannot constrain /^offset? the shear-only fit is not shown in that 
panel. Bins for which the x 2 value indicates a poor fit (> 30 for shear- and magnification-only fits, > 50 for the fits to both) are not 
shown, as the best-fit parameters for these bins are not meaningful. 





Lens Redshift Bin Centre 

Figure 7. The ratios of errors in our model parameters to the errors obtained from using shear information alone, for various lens 
mass and redshift bins. From left to right, the columns show the parameters: One-halo mass, Group mass, Satellite fraction, and the 
square-root of the determinant of the covariance matrix, which is representative of the total volume of the confidence regions. The red 
open circles show the ratios of errors from magnification alone to the errors from shear alone, the black stars show the ratios of errors 
from the combination of shear and magnification to shear alone, and the magenta crosses show the ratio of the theoretical combination 
of errors from magnification and shear to the errors from shear alone. Note that due to the multi-parameter nature of the fit and the fact 
that the best-fit values from shear and magnification alone sometimes disagree, the errors from the combined dataset are occasionally 
greater than for shear alone. 
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rameters at low redshift. However, as expected it is more 
beneficial at high redshift, particularly for the parameters 
relating to the host groups of the galaxies in the sample, 
M gr and / sa t. Interestingly, the errors of the fit to the com¬ 
bined dataset do not generally follow the theoretical calcu¬ 
lations, in some cases being even smaller, and in some cases 
being even larger than in the shear-only case. This suggests 
that the situation is more complicated than the theoretical 
calculation assumes. 

What is likely going on is a combination of two fac¬ 
tors. The shear and magnification signals and the fits to 
them have different degeneracies, and the two fits disagree 
on the best-fit set of parameters in some scenarios. The for¬ 
mer factor can result in the fit to the combined dataset hav¬ 
ing significantly smaller errors than the fit to either alone, as 
the combination of the two datasets breaks the degeneracies 
that are present in either alone, allowing much greater preci¬ 
sion. The latter factor can result in the fit to the combined 
dataset showing larger errors as the different regions pre¬ 
ferred by each dataset will spread out the most-reasonable 
region in the fit to the combined dataset. The result is that 
when the regions preferred by the individual fits overlap, 
the fit to the combined dataset is much smaller, while when 
they do not overlap, it can be larger. 

Notably, at high redshifts and low lens stellar masses, 
the shear signal approaches zero and isn’t reliably de¬ 
tectable, as can be seen in the bottom panel of Fig. [3] This 
results in the best-fit model parameters being unreliable 
when shear alone is used. However, as can be seen in the 
bottom panel of Fig. [S] the magnification signal can in fact 
still be measured for these bins, and so the best-fit model 
parameters from using only magnification are likely to be 
more reliable. 

In order to help illustrate the possible benefits of magni¬ 
fication data, Fig. [8] shows the confidence regions for the fit¬ 
ted parameters to two lens bins, selected to illustrate cases in 
which magnification is particularly useful. These show that 
magnification provides the most benefit in the M gr — / sa t 
parameter plane. Here, the shear signal is only able to pro¬ 
vide an upper boundary to the combination of M gr and / sa t. 
The magnification signal on its own is not much better, also 
covering larger regions of parameter space. However, the in¬ 
tersections of these confidence regions are relatively small, 
and so the fit to the combined dataset has a much tighter 
confidence region. 

It remains for future investigation to determine how 
much benefit magnification data will provide when all rel¬ 
evant biases are taken into account. The fact that even for 
this preliminary analysis, magnification occasionally pro¬ 
vides large benefits in reducing the errors in fitted param¬ 
eters is strongly encouraging and motivates further investi¬ 
gation to determine whether this will remain the case in a 
proper treatment. 


5 DISCUSSION 

In Section I5T1 we compare our method to the optimally- 
weighted correlation function method, which has been used 
in past analyses. In Section [5.21 we present a proposed im¬ 
plementation of our methodology for mass-mapping. In Sec¬ 
tion [El we present an estimate of the possible impact of 


source sample contamination on our results. In Section [El 
we test and discuss the impact of our choice to weight our 
measurement of the unmagnified source density on a per- 
field basis. Finally, in Section 15.51 we discuss further possible 
work. 


5.1 Comparison to Other Methods 


When applied to galaxy-galaxy lensing, the results of 
this methodology are nearly equivalent to those ob¬ 
tained through an optimally-weighted cross-correlation 
fu nction estimator. Th is method was first proposed by 
iMenard fe Bart elmannl (120021') , and in typical implementa¬ 
tions (eg. Scra nton et al. | 20051 : IHildebrandt et all l2009i ; 
iFord et al.ll2012i )~ uses a cross-correlation function in which 
all observed sources are weighted by a(m) — 1. The measured 
correlation function can then be compared to the predicted 
correlation functions from mass models. To date, it has not 
been stated in the literature how this correlation function 
relates to the magnification, and so we investigate this here. 

The angular cross-correlation represents the excess 
number density of galaxies in one sample in angular bins 
around galaxies in another sample. For scenarios in which 
only magnification is expected to influence the correlation 
between lens and source positions, which generally requires 
the two samples to be well-separated in redshift, the ex¬ 
pected cross-correlation between lenses and sources in a 
given angular and magnitude bin would be 


(■ w(0,m )) = 


No{m) 


- 1 = A t ( 60“ (m) - 1 - 1 , 


(39) 


where No(m) is the unmagnified number count of source 
galaxies at a given magnitude and ( N(9, m)) is the expected 
magnified number counts of galaxies at a given magnitude 
in a given angular bin relative to the lens sample, and us¬ 
ing Equation © to identify this ratio as corresponding to 

/x( 60 Q(m) -i. 

Approaches by eg. IScranton et al.l (l2005lf . 
IHildebrandt et all (120091 ') . and I Ford et alJ (120121 ') cal¬ 
culate the following optimally-weighted cross-correlation 
functi on, using a weighted version of the lLandv fe Szalavl 
rtl993f) estimator, 


opt 

U LS — 


LS“ -1 - RS° 


- (a - 1)LR 


RR 


+ <«-!>, (40) 


where S stands for sources, L stands for lenses, and R stands 
for the random catalogue, which is at least 10 times as large 
as the lens and source catalogues. All sources are weighted 
by their a(m) — 1 value, (a — 1) is the mean value of a(m) — 1 
across the entire sample, and all pair counts are normalized 
by the product of the total numbers of galaxies in each cat¬ 
alogue. 

This can be more easily compared to our method by 
instead presenting the cross-correlation function using the 
iDavis fe Peebles! j l983|) estimator, which converges to the 
same result as the lLandv fe Szalavl (Il993l f estimator at high 
number counts. The optimally-weighted cross-correlation 
function with this estimator is: 

t ™-i 

»D P r = %- -(«-!>, (41) 

using the same notation as in Equation ©3). The optimally- 
weighted correlation can therefore be seen as a measurement 
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Figure 8. Confidence regions of fitted parameters for shear-only, magnification-only, and combined fits to the 0.2 < z < 0.3, 1O 1O M0 < 
Aliens < 10 11 Mq bin (left) and the 1 < z < 1.1, 1O 1O M0 < M\ ens < 10 11 M@ bin (right). The thick and thin contours show 95% and 
68% confidence regions respectively. Note that the shear-only fits aren’t shown for the panels which include the /^offset parameter, as this 
parameter is unconstrained by the shear signal. 


of 

. opt £,ea„ n u H m j) - !] £,esamp H m i) - 1] (10) 

N 0 

,annu No ,samp 

where j £ samp represents the indices of all source galaxies 
in the sample, j £ annu represents only those source galax¬ 
ies which lie in the appropriate annulus behind any lens in 
the sample, Ao.samp is the total number of source galaxies 
in the sample, and iVo.annu is the expected number of source 
galaxies lying in the annuli in the presence of no magnifi¬ 
cation. If we convert this sum to an integral in the same 
fashion as in Equation m, we get: 

^o P t = fa n ( m ) \ a ( m ) - 1] dm - I* n 0 (m) H m ) - 1] dm 
f b no(m)dm 

(43) 

Compare this to the expansion of Equation 

n(m) [a(m) — 1] dm — f b no{m) [a(m) — 1] dm 
ff no(m) [a(m) — l] 2 dm 

(44) 



In both equations, the number count of galaxies ob¬ 
served within an annulus of interest only appears within 
one term, J n(m) \a(m,j) — 1] dm. The methods thus have 
equivalent statistical power, and it is possible to convert be¬ 
tween their measurements through: 

w opt 

A = 1 + 77-7T2T ’ ( 45 ) 

([«-!] ) 

and 


w opt = ([a- 1] 2 )(A- 1), (46) 


where 

, „ f b no(m) \a(m) — ll 2 dm 

([ Q ~ l] 2 ) — ] -• (47) 

J a no(m)am 

This comparison allows the statistical error of w opt to 
be easily estimated without expending computational power 
on a bootstrap or jackknife approach. It can instead be es¬ 
timated as: 


(J w opt 


No -d sairi r) 


(48) 


The derivation of this equation is presented in Equa¬ 
tion (IA2I) . This estimate relies on the assumption that back¬ 
ground galaxies follow an ideal Poisson distribution. Since 
clustering is expected among these galaxies, their actual dis¬ 
tribution will have a larger standard deviation than that of 
a Poisson distribution, and so this estimate will be biased 
low. 

Computationally, both our method and the correlation 
function method scale as at best 0(N\ log N s ) with the num¬ 
ber of lens galaxies N\ and number of source galaxies N s . 
This is due to the fact that they both include a step in which 
the source galaxies residing in a given annulus around lens 
galaxies must be determined, a process which can have this 
scaling at best. The difference comes in the coefficient of 
this term, which is dominated by the calculation of errors. 
Our method can calculate errors with only a few extra cal¬ 
culations per lens, as detailed in Section m However, for 
the correlation function method, a bootstrap or jackknife 
method is typically used to determine the errors, which re¬ 
quires re-running the analysis on s ubsets of the data multiple 
times (eg. 50 times in the case of lFord et alj|20l3 l. 
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The comparison we show here allows calculations of the 
magnification signal through a correlation function estima¬ 
tor to use Equation (0HJ as an estimate of the error. How¬ 
ever, this estimate will be biased low due to the clustering 
of source galaxies. The algorithm we use to calculate a more 
accurate error estimate for our method (see Equations 1281 
through 021) cannot be readily translated to a correlation 
function estimator, as it relies on the unmasked fractions of 
the annuli surrounding each lens in the sample, which are not 
calculated through a correlation-function-based approach. 
This results in our method being significantly faster than 
a correlation-function-based approach when accurate errors 
on the magnification estimate must be calculated, thanks to 
the fact that it does not have to resort to a computationally- 
expensive bootstrap or jackknife approach. 

It’s worth noting that many correlation function esti¬ 
mators, when measuring the function for large separations, 
bin objects i nto cells and determin e the unmasked fraction of 
each cell feg. lScranton et 21120051 ). The correlation function 
between these cells is then calculated through a summation 
over pairs of cells. In this case, the error can be determined 
empirically through an analogous calculation to the one used 
for our method. However, such binning is only practical at 
large separations. Modifications can be made to the small- 
scale estimation to allow for error calculations, but this will 
require the calculation of the unmasked fraction of each an¬ 
nulus. If this is done, it is then no longer necessary to use a 
random catalogue to determine the overdensity within each 
annulus; since its area is known, its overdensity can be more 
efficiently calculated by simply comparing its density to the 
density of the full observation. At this point, the method is 
in fact more similar to our proposed method than to the tra¬ 
ditional method of estimating a correlation function, lacking 
only the generalization done in our method to allow differ¬ 
ent redshifts for lenses. It is thus apparent that our proposed 
method is the most practical means to accurately determine 
the error in a measurement of a magnification signal. 

Finally, our method also provides the advantage 
that it can handle overlapping lens and source samples. 
Correlation-function-based approaches require the lens and 
source samples to be well-separated in redshift space so that 
there is no expected clustering between the two samples. Our 
method has the benefit that each lens is treated individually, 
and so it is only necessary that each lens be well-separated 
from the subset of the source galaxy sample which we use to 
measure magnification. This results in the ability for lenses 
at the low-redshift end of the lens sample to use a greater 
amount of redshift space to measure the magnification sig¬ 
nal, regardless of the fact that higher-redshift lenses in the 
sample cannot use this space. This will result in an increased 
signal-to-noise of the magnification measurement. 

5.2 Proposed Implementation for Mass Mapping 

As we showed in Section l5.ll our method has equivalent sta¬ 
tistical power to cross-correlation function methods already 
in use when applied to galaxy-galaxy lensing. One important 
advantage this method has is that it is capable of provid¬ 
ing an estimate of the mean magnification of an arbitrary 
region of space, not limited to a circular annulus, or even to 
a region of constant weight. 

As the statistical power of magnification varies with 


the square-root of the area of the region measured, mass 
mapping with magnification alone must necessarily involve a 
trade-off between resolution and accuracy. The simplest ap¬ 
proach to mass mapping would involve measuring the mag¬ 
nification, and thus projected overdensity, in a circular aper¬ 
ture around any point in the region of interest. While this 
method can be easily applied with the methodology outlined 
in the previous sections, it suffers from the drawback that 
the resulting mass map will show discontinuous behaviour 
whenever a source galaxy enters or exits the aperture. 

This discontinuous behaviour can be avoided through 
instead applying a filter function which approaches zero at 
large separation, such as a two-dimensional Gaussian weight 
function. This also provides the benefit that source galaxies 
closer to the point of interest will have more influence on 
its measured signal than distant galaxies. Let us call this 
function F(9). 

The magnification at a given angular position 8 and 
redshift £ can then be estimated as 

Uki 1 + + ^ 2 ) x 

^F — 9^j [a(rrij) — 1] — Apd(z + , 

(49) 

where d(z) and w(z) are defined as in Section [2.21 j repre¬ 
sents a summation over all source galaxies at redshift z + Az 
or greater, with 8j representing the angular coordinates of 
these galaxies, and 

Ap = Jj F(8)dSl. (50) 

In practice, source galaxies sufficiently distant from the po¬ 
sition of interest such that F — (?j j —» 0 can be ignored 
in this calculation. 

One notable issue with mass-mapping with magnifica¬ 
tion is the fact that position-dependent biases in the redshift 
determinations of source galaxies can result in a spurious 
magnification signal, as discussed in Section \2 . 3 1 In the sce¬ 
nario of galaxy-galaxy lensing, this effect is of large enough 
scale that it can be treated as a flat offset and included 
as a fitted parameter in the model, but this is not possible 
for mass mapping. The effect is somewhat reduced in mass¬ 
mapping, as the positions of lenses are not used, and so this 
effect should only appear in tomographic analyses in which a 
redshift cut-off is imposed on sources. This cut-off would al¬ 
low for position-dependent effects, as the number of sources 
beyond the cut-off will depend on the errors in their redshift 
determinations, which is itself position-dependent. Care will 
thus have to be taken in such an analysis to account for this 
effect. 

5.3 Estimating Source Sample Contamination 

As shown in Fig. [3] the fraction of galaxies whose photo¬ 
metric redshifts make them appear to be sources but which 
in fact lie within the lens plane can is of order 0.1 per cent 
at low redshifts, and can be as high as ~ 10 per cent at 
high lens redshifts, even with a large buffer redshift between 
lenses and sources imposed. Galaxies in the lens plane will 
typically be correlated in position with each other due to 
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Lens Redshift Bin Centre 

Figure 10. The change in the fitted values of re 0 fj- aet when weight¬ 
ing is applied on a per-field basis in the source-counting proce¬ 
dure. This weighting accounts in part for observational effects 
such as variable seeing which indirectly result in the number den¬ 
sities of lenses and sources being correlated or anti-correlated. The 
magnitude of the change shown here provides an estimate of the 
uncertainty in the value of re 0 ff set due to smaller-scale variations. 

their mutual gravitational attraction. Therefore, if galaxies 
in this plane are misidentified as sources, we would expect 
them to generally appear closer to lens galaxies than they 
would through random chance alone. 

We can make an order-of-magnitude estimate of the im¬ 
pact this will have on our measured magnification signal 
through multiplying the correlation function of galaxies in 
the lens plane with the contamination and a factor repre¬ 
senting the conversion from the correlation function to re. 
For this conversion factor, we choose 0.25, which would cor¬ 
respond to a mean | a(m) — 1| of ~ 0.5. Typically, the mean 
is less than this, so this results in an upper bound on the 
contamination of re. We plot the resulting predicted con¬ 
tamination for selected lens redshift and stellar mass bins in 
Fig. 0 As can be seen through comparison with Fig. [oj this 
value is typically negligible for z < 0.6, becoming compara¬ 
ble in magnitude to re in the 0.6 < z < 0.7 bin. In higher 
redshift bins, the effect is even worse. 

Further investigation of correcting for this effect is be¬ 
yond the scope of this paper, but it is unlikely to be present 
significant difficulties. By virtue of the fact that the galaxies 
in the lens plane which are assigned photometric redshifts 
appropriate for sources, they likely have a similar magnitude 
distribution to legitimate source galaxies, and so correction 
for their differing magnitude distribution may not be neces¬ 
sary. As such, a simple calculation based on the correlation 
function of galaxies in the lens plane and the mean value of 
a — 1 for sources may be all that is needed. 

5.4 Effects of Field Weighting 

In Section we discussed how correlations and anti¬ 
correlations of lenses and sources in regions of the sky can 
result in an effect which mimics a magnification signal. 
We chose to partially correct for this effect by applying 
weights on a field-by-field basis in determining the unmagni¬ 


fied source density. We can demonstrate the impact of this 
weighting through a comparison of the best-fit re 0 ffset pa¬ 
rameters in the scenarios in which we do and do not apply 
field-by-field weighting, as shown in Fig. 1101 

As expected, the application of weighting has the same 
effect on all lens stellar mass bins, but the effect varies with 
redshift. This is consistent with our hypothesis that variable 
seeing could result in varying amounts of scatter between 
redshift bins from field to field, as this effect would not be 
expected to be constant across all redshifts due to the com¬ 
plicated nature of determining photometric redshifts. 

In general, the magnitude of change in re 0 ff se t due to 
field-by-field weighting is < 0.004. Any remaining effect due 
to variable seeing on smaller scales is likely to result in a 
smaller change, and so we would expect the proper fitted 
values for re 0 ff se t for all bins to be of smaller magnitude. From 
the fitted values shown in Fig. [6] it appears that the scatter 
in fitted values of re 0 ff se t is comparable to this magnitude. 
Additionally, many bins of the same redshift but different 
stellar masses show significantly different fitting values of 
^offset- This suggests that the fitted values for re 0 fF set are not 
sensitive to this effect alone, but are also affected by other 
aspects of the measured lensing signals. This is possibly due 
to the fact that the model we use does not include contri¬ 
butions for larger-scale structure, which is dominant at the 
same large radial separations that re 0 ff se t is. 

Alternatively, let us consider that this may be due not 
to an observational effect, but due to actual over- or un¬ 
derdensities in the regions of these lense s. The convergence 
maps generated by [Van Waerbeke et al.1 (i2013i ) from shear 
data alone show that, when smoothed on scales of 2' and 
correcting for the contribution of noise, the convergence has 
typical standard deviation of ~ 0.0060 The typical area cov¬ 
ered by all lenses within one of our stellar mass and redshift 
bins is of order ~ 10 5 arcmin 2 . As this is roughly four orders 
of ma gnitude greate r than the smoothing scale applied by 
IVan Waerbeke et al.l (120131 ). we would expect the standard 
deviation of re. due to large scale structure within this region 
to drop by two orders of magnitude, to ~ 0.00006. Clearly, 
this is too small to account for the scatter we observe in the 
fitted values of re 0 ff se t, and so they are more likely due to 
observational effects or fitting errors than cosmic variance. 

5.5 Further Research 

In our attempts to fit a mass profile to shear and magnifica¬ 
tion data simultaneously, we found several indications that 
our fit was not adequate, including poor \ 2 comparisons, 
extreme and unlikely values fit for certain bins, and incon¬ 
sistent trends over changing lens redshift and stellar mass. 
It is possible that these issues might be resolved through the 
use of a different model mass profile, or through the inclusion 
of the effects of large-scale structure on the lensing signal. 

6 We determi ned this value through mas s map data provided by 
the authors of lVan Waerbeke et al.l ll2013l l. Under the assumption 
that the variance in the predicted convergence from their mass 
maps is the sum of the variances due to noise and from the actual 
variance in re, we took the difference between the variances in re 
for their data maps and noise maps to estimate the variance, and 
thus standard deviation, in re due to the actual distribution of 
matter. 
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Figure 9. The maximum change in the measured k due to contamination of the source sample with lenses. Low-redshift bins are shown 
in the top panel, and high-redshift bins are shown in the bottom panel. This maximum change is calculated from the product of the 
correlation function shown in Fig. EH the contamination at the lens redshift, which is shown in Fig. [3] and the factor 0.25, which 
represents a conservative upper bound on the slope of k/w. 


This will allow us to better determine the potential benefits 
of magnification in breaking degeneracies between parame¬ 
ters characterizing the environment surrounding galaxies in 
our samples. 

As mentioned in Section m dust obscuration can re¬ 
sult in the measured magnification signal to be spuriously 
high if it is not corrected for. In this paper, we handled 
this effect through discarding the two innermost radial bins 
when we fit our mass profiles to the magnification data, as 


these are the bins which are likely to show the most signif¬ 
icant contribution from dust. A proper analysis would re- 
quire modellin g the obscuration from dust, as done by eg. 
iMenard et al.l ( 2010T1 . This will be necessary before firm con¬ 
clusions can be drawn from our model fits. 
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6 SUMMARY AND CONCLUSIONS 

In this paper, we have presented a new method to estimate 
weak lensing magnification which can be applied to an ar¬ 
bitrary region of space, as expressed in Equation ca and 
in a more efficient form for computations in Equation (1231) . 
We’ve shown that the estimator resulting from this method 
has equivalent statistical power to the optimally-weighted 
correlation function estimator which has been used for past 
analyses, and we provide conversions between them in Equa¬ 
tion (H51) and Equation m- This conversion also allows a 
simple way to estimate the errors on measurements made 
with the correlation function estimator without the use of 
bootstrap or jackknife methods, which we present in Equa¬ 
tion on. However, this estimate is likely to be biased low 
due to the clustering of source galaxies. Our method can 
calculate a more accurate estimate of error with only a few 
extra calculations, which presents a large benefit over the 
requirement of correlation-function-based methods to use 
bootstrap or jackknife approaches to accurately calculate er¬ 
ror. Our method is also capable of handling overlapping lens 
and source samples, which can increase the signal-to-noise 
of the magnification signal when this is the case. 

The estimator we have presented here provides a suit¬ 
able mathematical basis to use magnification for mass map¬ 
ping, for instance through the form in Equation (1491) . How¬ 
ever, applying this will be difficult in practice, as the con¬ 
tributions to the magnification signal from actual fluctua¬ 
tions in the projected mass density field are degenerate with 
localized observational biases, such as due to variable see¬ 
ing across the field of view. It remains for future research 
to determine whether this issue can be resolved, through 
advances in photometric redsliift estimation to correct for 
seeing-induced biases or other techniques. 

We’ve shown a proof-of-principle of our method applied 
to galaxy-galaxy lensing in the CFHTLenS fields, and shown 
that the convergence profiles determined with it are con¬ 
sistent with the convergent predictions from mass profiles 
which are fit to shear data in the fields, as can be seen in 
Fig. [5] As expected, in most cases magnification data pro¬ 
vides poorer signal-to-noise than does shear, but its con¬ 
tribution is not always negligible. In fact, in certain cases 
its combination with shear seems to provide significantly 
smaller errors on model parameters than either alone, likely 
due to the magnification data helping break degeneracies 
between multiple parameters. This finding is tentative, how¬ 
ever, as our present algorithm was only designed as a proof- 
of-principle and does not yet include provisions for effects 
such as dust obscuration, which can result in biases in the 
magnification signal. 
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APPENDIX A: DERIVATIONS OF SELECTED 
EQUATIONS 


In Equation 0, we presented the expected standard devia¬ 
tion of the estimator /r i for the magnification of an individ¬ 
ual magnitude bin. This equation can be derived as follows, 
starting from Equation 0: 


O dni 

| OLi — 1| dn 0 ,i 

y/dn,i 

|Oi — 1| dn 0 ,i 
Cdno.ifJ . 0 - 1 
\cti — 1| dno,i 
\/l + — l)(/r — 1) 

| on - 1| y/dno,i 


(Al) 


In Equation (1481) , we presented an equation for estimat¬ 
ing the standard error of the optimally-weighted correlation 
function estimator for magnification, without relying on a 
bootstrap or jackknife approach. This equation can be de¬ 
rived as follows, starting from Equation (|46j: 


<j w opt — ([a — 1] )o£ 


no (to) [ a(m) — l ] 2 dm 
no(m)dm 

x (^J n o (to) [ ck ( to ) — l ] 2 dm 

( la n ° (to) [a{m) — l ] 2 dm j 
no(m)dm 


- 1/2 


1/2 


<[«-l] 2 > 


No -Asamp 

where No is defined as in Section I5T1 


(A2) 


APPENDIX B: SUPPLEMENTAL FIGURES 

Fig. [BT] shows the correlation function between galaxies in 
the lens and source samples within the same redshift slices 
for low-redshift bins. This can be used to estimate the con¬ 
tamination fraction of lenses in the source sample, which we 
showed in Fig. [3] 


APPENDIX C: CODE ACCESS 

In order to facilitate further work with magnification, 
we make the magnification and shear measurement 
code used for this paper publicly-available. The code 


© 2015 RAS, MNRAS 000. HIM 













20 B.R. Gillis and A.N. Taylor 



Projected Separation (kpc) 



0 500 1000 1500 2000 500 1000 1500 2000 500 1000 1500 2000 500 1000 1500 2000 500 1000 1500 2000 500 1000 1500 2000 


Projected Separation (kpc) 


Figure Bl. Correlation function between galaxies from the lens and source samples in various redshift and stellar mass bins. Low 
redshift bins are shown in the top panel, and high redshift bins in the bottom panel. 


may be accessed online through a repository hosted at 
https://bitbucket.org/brgillis/magnification_public. 
If this repository is no longer available and a web search 
fails to locate a new host, interested parties are invited to 
contact the lead author, Bryan Gillis, at brg@roe.ac.uk to 
request the code. 
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